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Abstract 

The original Lindhard-Scharff-Schi0tt (LSS) theory and the more recent Tilinin theory for calculating the nuclear 
and electronic stopping powers of slow heavy ions are compared with predictions from the SRIM code by Ziegler. 
While little discrepancies are present for the nuclear contribution to the energy loss, large differences are found in 
the electronic one. When full ion recoil cascade simulations are tested against the elastic neutron scattering data 
available in the literature, it can be concluded that the LSS theory is the more accurate. 
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1. Introduction 

It is well known that an ion moving inside a 
medium can loose energy by collisions with both 
electrons and nuclei. Theoretically, it is necessary 
to account for both processes to reach an accu- 
rate description of the energy loss below a few 
keV/amu. From the experimental point of view, 
most detectors are just sensitive to electronic en- 
ergy loss. Knowledge of the energy loss sharing 
between the two processes is mandatory for a de- 
tailed understanding of the response to particles 
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interacting with the detecting medium through nu- 
clear recoils. This is the case for Weakly Interact- 
ing Massive Particles (WIMPs) , which are possible 
constituents of the galactic dark matter. 

The purpose of the present work is to explore the 
different available descriptions for the two parts 
of the energy loss between 1 and 100 keV. Range 
measurements are difficult at such low energies and 
almost all published data have been obtained em- 
ploying elastic neutron scattering for transferring 
small and known amounts of energy to atoms of the 
detecting material. To reduce the theoretical diffi- 
culties, the attention is focused on pure substances: 
hence only symmetric projectile/target atom com- 
binations will be investigated. From all the avail- 
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able measurements, known to the authors, the fol- 
lowing are then selected: Si [3,4], Ge [5,6], and liq- 
uid Xe [7,8,9,10]. Liquid Ar will also be considered 
on account of its interest for dark matter searches. 

Among the few calculations from first principles 
without free parameters, there are still the origi- 
nal theory of Lindhard [1] and its reevaluation by 
Tilinin [2]; they will be examined here. 

From all the available codes, only SRIM [11] will 
be discussed, because: i) it gives separately nuclear 
and electronic energy losses, ii) it covers the low 
energy range of interest, iii) from an independent 
survey [12] it was found the most accurate. 




Ion energy E [keV] 

Fig. 1. Nuclear stopping power Sn as a function of the ion 
energy for the symmetric projectile/target combinations 
considered in the present study. 



2. The nuclear stopping power 

According to the Lindhard-ScharfT-Schi0tt 
(LSS) theory [1], the nuclear stopping power S„ 
of a heavy ion is best described by rescaling its 
energy E and range R to the non-dimensional 
variables e and p, respectively, defined as 
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with Z = Z 2 p /3 + Z T '° and A tot = Ap + At- In 
Eq. (1) N is the number density of the target ma- 
terial, Zp and Ap are the atomic and atomic mass 
numbers of the projectile atom, respectively, Zt 
and At are the correspondents for the target atom, 
<zb is the Bohr radius, E-& the Bohr energy and Ctf 
the Thomas-Fermi constant (97r 2 /2 7 ) 1 / 3 . In fact, 
while the nuclear part of the energy loss (dE/dx) n 
depends on the projectile and target ions, (de/dp) n 
depends on the Thomas-Fermi interaction poten- 
tial alone and is a universal function /(e), which 
can be calculated numerically [13]. The values of 
S n are reproduced with the physical units in Fig. 1. 

A similar rescaling is employed in SRIM (Z in 
all appearances in Eq. (1) is replaced by an empir- 
ical Z 1 ! 2 = Z p 23 + Z T - 23 ), but now the universal 
function /(e) is determined not from first princi- 
ples but from a fit to experimental data. The re- 
sults are also shown in Fig. 1. 
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While differences are small, it is known from 
sputtering data [14] that the LSS theory overesti- 
mates S n . This problem is reduced in SRIM, as it 
is based on data themselves. Lindhard et al. [13] 
had also warned that the Thomas- Fermi treatment 
might be inadequate for very low energies (i.e. e < 
10 -2 ), where mostly the tails of the ion-ion poten- 
tial are probed. In the present study, this is of con- 
cern only for Xe nuclear recoils below 15 keV. 



3. The electronic stopping power 

The electronic energy loss (dE/dx) e of a unitary 
charge particle with a velocity (3 was described 
as an interaction with an electron plasma in the 
original work of Fermi and Teller. They explicitly 
distinguished two cases, for /3 above and below 
the Fermi velocity (3p. For f3 < /3p, (dE/dx) e was 
found to be proportional to /?, with a proportion- 
ality coefficient being a unique function of the elec- 
tron plasma density no , usually expressed in terms 
of the Wigncr- Seitz radius r 3 = (3n /(4 7 r)) 1 / 3 [15]. 
Their result is plotted in Fig. 2. Typically, r s needs 
to be corrected because tightly bound electrons 
contribute only marginally to (dE/dx) e much be- 
low the Bragg peak. In crystals, like Si and Ge, this 
effective density of the free electron plasma can be 
deduced from optical properties [16]. In liquids, 
like Ar and Xe, the problem is much more difficult 
and no correction was attempted in Fig. 2. Succes- 
sively, Lindhard calculated in a self consistent way 
the local increase of the electron plasma density 
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around the intruder particle due to its Coulomb 
field. This leads to a higher (dE/dx) e as can be 
seen in Fig. 2. Finally, Ritchie considered the case 
where the Coulomb field is exponentially screened, 
slightly decreasing (dE/dx) e (see Fig. 2). The last 
effect is of particular relevance for ions, which 
can accommodate bound states while sweeping 
through the electron plasma. For a bare ion, a 
scaling with Z P to the elementary particle case is 
expected, as assumed in Fig. 2. In reality, this is 
not correct and the theory was extended to a par- 
tially ionized intruder by Fcrrcll and Ritchie [15], 
but the determination of the equilibrium charge 
of a given ion remains a difficult task. Lindhard 
also independently investigated this problem [1] 
and, using the Thomas-Fermi theory, arrived at 
a closed form for the proportionality coefficient. 
In terms of the non-dimensional variables intro- 
duced in Eq. (1), his result can be expressed as 
(de/dp) e = k \fe where 
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with £ e m z\J^ (which is regarded only as an ap- 
proximation by Lindhard) [1]. The points corre- 
sponding to the projectile/target combinations of 
interest for the present study are also reported in 
Fig. 2. The suppression of the electronic energy loss 
occurs mostly due to the partial ionization of the 
intruder. It strongly increases with Zp (Zp = Zp). 

The proportionality of (dE/dx) e with (3 is also 
a feature of the SRIM code, allowing a value of 
(dE/dx) e /(/3 Z P ) to be extracted (see Fig. 2). 
While SRIM exceeds the LSS theory for Si, it then 
decreases consistently below it, up to a factor of 
4 for Xe. The reason for this discrepancy is un- 
clear, because the details on the implementation 
of (dE/dx) e in SRIM for low velocities are not 
public, but it probably resides in the estimate of 
the intruder charge state. It has been verified that 
data for protons of comparable energies per nu- 
cleoli on Ar and Xe are well reproduced by SRIM. 

The big drawback of the described theoretical 
approaches is to assume that the electronic and 
nuclear collisions are uncorrelated [1]. In reality, 
the screened Coulomb repulsion between the two 
interacting nuclei makes part of the impact pa- 
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Fig. 2. Proportionality coefficient of the electronic energy 
loss to the particle velocity as a function of the Wigner- 
Seitz radius r s . For ions a Zj, scaling is applied. 



rameter range unavailable for the scattering of 
the electrons belonging to the target atom in the 
screened Coulomb field of the projectile and vice 
versa. Tilinin [2] has shown that the final net effect 
is a great decrease of S e for e <C 1, with a corre- 
sponding lack of proportionality to -y/e. His results 
can also be recast in the form of Eq. (2) where £ e 
is replaced by a function t(c, Zp/Zp) that can be 
tabulated [2]. In the present case, his theory pre- 
dicts roughly half the value of S e expected from 
LSS and SRIM for Si and approximately agrees 
with SRIM for Xe. 



4. The full ion recoil cascade 

As mentioned, resort must be made to response 
measurements employing elastic neutron scatter- 
ing, where recombination or quenching may influ- 
ence the final fraction of the total energy trans- 
ferred to electrons that is detectable as excitation 
or ionization, particularly in the case of scintilla- 
tion yield for LXc. Theoretically, however, an even 
bigger disadvantage is present; especially for high 
Z elements: S n dominates over S e and most of the 
primary ion energy is transferred to nuclear re- 
coils. The knocked ion undergoes the same process, 
resulting in a full cascade of recoils, whose total 
electronic energy loss must be evaluated. In the 
case of the SRIM code, a second program, called 
TRIM [11], reads S n and S e from the first and per- 
forms the computation. Lindhard and his group 
solved numerically the transport equations corre- 
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Fig. 3. Total fraction of the initial ion energy transferred 
to electrons integrated over the full cascade as a function 
of the ion energy itself. 

sponding to the LSS theory and arrived at a pa- 
rameterization in terms of e and n (see Eq. (9) of 
Ref. [17]). Both results are compared with data 
in Fig. 3 (the two series of points for Ge are in- 
consistent). The LSS theory is on average better, 
while SRIM both ovcrpredicts and strongly undcr- 
predicts the data for Si and Xe, respectively. While 
the first discrepancy could be ascribed to recom- 
bination, the second appears even more surprising 
due to the probable presence of quenching. 



5. Conclusions 

For slow heavy ions, the nuclear stopping power 
predicted by the original LSS theory and the cur- 
rent SRIM code differ at most by « 15%. On the 
contrary, for the electronic stopping power, big dis- 
crepancies are present between the LSS theory, the 
theory of Tilinin and SRIM (up to a factor of sa 4 
for Xe) . Judging from the elastic neutron scatter- 
ing data, the LSS theory seems the best of all. More 
detailed full cascade simulations will be performed 
in the future to assess the robustness of this conclu- 
sion. New experimental data for Germanium would 
be highly needed for clarifying the situation. 
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